Quantum Simulation of Molecules in Solution.

Quantum chemical calculations on quantum computers have been focused mostly on simulating molecules in the gas phase. Molecules in liquid solution are, however, most relevant for chemistry. Continuum solvation models represent a good compromise between computational affordability and accuracy in describing solvation effects within a quantum chemical description of solute molecules. In this work, we extend the variational quantum eigensolver to simulate solvated systems using the polarizable continuum model. To account for the state dependent solute-solvent interaction we generalize the variational quantum eigensolver algorithm to treat non-linear molecular Hamiltonians. We show that including solvation effects does not impact the algorithmic efficiency. Numerical results of noiseless simulations for molecular systems with up to 12 spin-orbitals (qubits) are presented. Furthermore, calculations performed on a simulated noisy quantum hardware (IBM Q, Mumbai) yield computed solvation free energies in fair agreement with the classical calculations.

Quantum Algorithm for Simulating Single-Molecule Electron Transport.

An accurate description of electron transport at a molecular level requires a precise treatment of quantum effects. These effects play a crucial role in determining the electron transport properties of single molecules, which can be challenging to simulate classically. Here we introduce a quantum algorithm to efficiently calculate electronic current through single-molecule junctions in the weak-coupling regime. We show that a quantum computer programmed to simulate vibronic transitions between different charge states of a molecule can be used to compute electron-transfer rates and electronic current. In the harmonic approximation, the algorithm can be implemented using Gaussian boson sampling devices, which are a near-term platform for photonic quantum computing. We apply the algorithm to simulate the current and conductance of a magnesium porphine molecule. The algorithm provides a means for better understanding the mechanism of electron transport at a molecular level, which paves the way for building practical molecular electronic devices.

Quantum algorithm for simulating molecular vibrational excitations.

The excitation of vibrational modes in molecules affects the outcome of chemical reactions, for example by providing molecules with sufficient energy to overcome activation barriers. In this work, we introduce a quantum algorithm for simulating molecular vibrational excitations during vibronic transitions. We discuss how a special-purpose quantum computer can be programmed with molecular data to optimize a vibronic process such that desired modes get excited during the transition. We investigate the effect of such excitations on selective bond dissociation in pyrrole and butane during photochemical and mechanochemical vibronic transitions. The results are discussed with respect to experimental observations and classical simulations. We also introduce quantum-inspired classical algorithms for simulating molecular vibrational excitations in special cases where only a limited number of modes are of interest.

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems.

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light-matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

Nonequilibrium Solvent Polarization Effects in Real-Time Electronic Dynamics of Solute Molecules Subject to Time-Dependent Electric Fields: A New Feature of the Polarizable Continuum Model.

We develop an extension of the time-dependent equation-of-motion formulation of the polarizable continuum model (EOM-TDPCM) to introduce nonequilibrium cavity field effects in quantum mechanical calculations of solvated molecules subject to time-dependent electric fields. This method has been implemented in Octopus, a state-of-the-art code for real-space, real-time time-dependent density functional theory (RT-TDDFT) calculations. To show the potential of our methodology, we perform EOM-TDPCM/RT-TDDFT calculations of trans-azobenzene in water and in other model solvents with shorter relaxation times. Our results for the optical absorption spectrum of trans-azobenzene show (i) that cavity field effects have a clear impact in the overall spectral shape and (ii) that an accurate description of the solute shape (as the one provided within PCM) is key to correctly account for cavity field effects.

Predicting signatures of anisotropic resonance energy transfer in dye-functionalized nanoparticles.

Resonance energy transfer (RET) is an inherently anisotropic process. Even the simplest, well-known Förster theory, based on the transition dipole-dipole coupling, implicitly incorporates the anisotropic character of RET. In this theoretical work, we study possible signatures of the fundamental anisotropic character of RET in hybrid nanomaterials composed of a semiconductor nanoparticle (NP) decorated with molecular dyes. In particular, by means of a realistic kinetic model, we show that the analysis of the dye photoluminescence difference for orthogonal input polarizations reveals the anisotropic character of the dye-NP RET which arises from the intrinsic anisotropy of the NP lattice. In a prototypical core/shell wurtzite CdSe/ZnS NP functionalized with cyanine dyes (Cy3B), this difference is predicted to be as large as 75% and it is strongly dependent in amplitude and sign on the dye-NP distance. We account for all the possible RET processes within the system, together with competing decay pathways in the separate segments. In addition, we show that the anisotropic signature of RET is persistent up to a large number of dyes per NP.

Modeling solvation effects in real-space and real-time within density functional approaches.

The Polarizable Continuum Model (PCM) can be used in conjunction with Density Functional Theory (DFT) and its time-dependent extension (TDDFT) to simulate the electronic and optical properties of molecules and nanoparticles immersed in a dielectric environment, typically liquid solvents. In this contribution, we develop a methodology to account for solvation effects in real-space (and real-time) (TD)DFT calculations. The boundary elements method is used to calculate the solvent reaction potential in terms of the apparent charges that spread over the van der Waals solute surface. In a real-space representation, this potential may exhibit a Coulomb singularity at grid points that are close to the cavity surface. We propose a simple approach to regularize such singularity by using a set of spherical Gaussian functions to distribute the apparent charges. We have implemented the proposed method in the Octopus code and present results for the solvation free energies and solvatochromic shifts for a representative set of organic molecules in water.

Scaling in the correlation energies of two-dimensional artificial atoms.

We find an unexpected scaling in the correlation energy of artificial atoms, i.e., harmonically confined two-dimensional quantum dots. The scaling relation is found through extensive numerical examinations including Hartree-Fock, variational quantum Monte Carlo, density functional, and full configuration interaction calculations. We show that the correlation energy, i.e., the true ground-state total energy minus the Hartree-Fock total energy, follows a simple function of the Coulomb energy, confinement strength and number of electrons. We find an analytic expression for this function, as well as for the correlation energy per particle and for the ratio between the correlation and total energies. Our tests for independent diffusion Monte Carlo and coupled-cluster results for quantum dots-including open-shell data-confirm the generality of the scaling obtained. As the scaling also applies well to ≳100 electrons, our results give interesting prospects for the development of correlation functionals within density functional theory.

Modeling opto-electronic properties of a dye molecule in proximity of a semiconductor nanoparticle.

A general methodology is presented to model the opto-electronic properties of a dye molecule in the presence of a semiconductor nanoparticle (NP), a model system for the architecture of dye-sensitized solar cells. The method is applied to the L0 organic dye solvated with acetonitrile in the neighborhood of a TiO2 NP. The total reaction potential due to the polarization of the solvent and the metal oxide is calculated by extending the polarizable continuum model integral equation formalism. The ground state energy is computed by using density functional theory (DFT) while the vertical electronic excitations are obtained by time-dependent DFT in a state-specific corrected linear response scheme. We calculate the excited state oxidation potential (ESOP) for the protonated and deprotonated forms of the L0 dye at different distances and configurations with respect to the NP surface. The stronger renormalizations of the ESOP values due to the presence of the TiO2 nanostructure are found for the protonated dye, reaching a maximum of about -0.15 eV. The role of protonation effect is discussed in terms of the atomic Löwdin charges of the oxidized and reduced species. On the other hand, we observed a weak effect on the L0 optical excitation gap due to the polarization response of the NP.

Selection of deep brain stimulation candidates in private neurology practices: referral may be simpler than a computerized triage system.

OBJECTIVE: The objective of this study is to compare a computerized deep brain stimulation (DBS) screening module (Comparing Private Practice vs. Academic Centers in Selection of DBS Candidates [COMPRESS], NeuroTrax Corp., Bellaire, TX, USA) with traditional triage by a movement disorders specialized neurologist as the gold standard. METHODS: The COMPRESS consists of a combination of the Florida Surgical Questionnaire for Parkinson disease (FLASQ-PD), a cognitive assessment battery provided by MindStreams® (NeuroTrax Corp.), and the Geriatric Depression Scale and the Zung Anxiety Self-Assessment Scale. COMPRESS resulted in the classification of patients into three categories: "optimal candidate,""probable candidate," and "not a good candidate." Similar categorical ratings made by a referring private practice neurologist and by a trained movement disorders specialist were compared with the ratings generated by COMPRESS. RESULTS: A total of 19 subjects with Parkinson's disease were enrolled from five private neurological practices. The clinical impressions of the private practice neurologist vs. those of the movement disorders specialist were in agreement approximately half the time (10/19 cases). The movement disorders specialist and COMPRESS agreed on 15/19 cases. A further comparison between outcomes from the entire COMPRESS module and the FLASQ-PD questionnaire by itself resulted in high agreement (18/19 cases in agreement). CONCLUSIONS: The COMPRESS agreed with an in-person evaluation by a movement disorders neurologist approximately 80% of the time. The computerized COMPRESS did not provide any screening advantage over the short FLASQ-PD paper questionnaire. Larger studies will be needed to assess the utility and cost effectiveness of this computerized triage method for DBS.